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Abstract 



We study the number of collisions X„ of an exchangeable coalescent with multi- 
ple collisions (A-coalescent) which starts with n particles and is driven by rates 
determined by a finite characteristic measure v{dx) = x~'^A{dx). Via a coupling 
technique we derive limiting laws of Xn , using previous results on regenerative com- 
positions derived from stick-breaking partitions of the unit interval. The possible 
limiting laws of Xn include normal, stable with index 1 < a < 2 and Mittag-Leffler 
distributions. The results apply, in particular, to the case when z/ is a beta(a — 2, b) 
distribution with parameters a > 2 and 6 > 0. The approach taken allows to de- 
rive asymptotics of three other functionals of the coalescent, the absorption time, 
^ the length of an external branch chosen at random from the n external branches, 

^vq ■ and the number of collision events that occur before the randomly selected external 

Tij" , branch coalesces with one of its neighbours. 

Keywords: absorption time, coupling, exchangeable coalescent, external branch length, 
(^ . multiple collisions, number of collisions, regenerative composition 
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1 Introduction 

5^ ■ Let n„ = (n„(t),t > 0) be an exchangeable coalescent with multiple collisions, as intro- 

duced by Pitman [22] and Sagitov [25]. This is a cadlag Markov process that starts at 
t = with n particles and evolves according to the following dynamics. When m particles 
are present, each fc-tuple of them collides and merges to form a single particle at rate 

\mk= f x''{l-xr-''u{dx), 2<k<m (1) 

J[o,i] 
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The definition of the rates requires that the characteristic measure v had a finite second 
moment. These processes are sometimes called 'lambda-coalescents', following Pitman's 
paper [22] where the rates were introduced by means of the finite measure A(dx) = 

The coalescent starting with n particles is usually understood as a process with values 
in the set of partitions of \n\ = {1, . . . ,n}. To that end, one needs to label the initial 
particles by integers and, tracking their evolution, to record which of the initial particles 
comprise a merger. In this context, the exchangeability means that the restriction of the 
process to [m] (or any other subset of [n] with m < n elements) is a stochastic copy of 
the process starting with m particles, see p] [23] for survey and details. For most of our 
purposes, however, it is sufficient to identify the coalescent with a pure death process of 
counting the particles, with state-space [n] and transition rates 

gn,m = ( , I \n,n-m+l, I < m < H - I 

\m — 1/ 
so that the total transition rate from state n is 

n— 1 „ 

9n = y^ gn,m = (1 - (1 - x)" - nx{l - x)"'^)z/(dx). 

m=l -^m 

In terms of the richer partition- valued process, the event n„(t) = m occurs when at time 
t the partition of [n] has m blocks, and a merge of some k > 2 blocks into a single block 
corresponds to a transition from m to m — k + 1. 

The coalescent starts in the state n„(0) = n and terminates in the state 1, as all 
particles eventually merge in one. The number of collisions X„ is the number of transitions 
which the process needs to proceed from n to the terminal state 1. The absorption time 
Tn = inf{t > : n„(t) = 1} is the time when the process enters the terminal state. Note 
that in biological applications of the coalescent theory the absorption time r„ has the 
meaning of time back to the most recent common ancestor. The main focus of this paper 
is on the weak asymptotic behaviour of X„ and r„ as n tends to infinity. In particular, 
the main result (Theorem 13.11) is a criterion describing all possible limit laws of, properly 
normalized and centered, X„: these can be normal, stable with index 1 < a < 2 and 
Mittag-Leffler distributions. Another result (Theorem 14. 2p provides sufficient conditions 
under which Tn possesses the same collection of limit laws. The limiting results for both, 
Xn and Tn, will be obtained by comparison with a simpler process which was partly 
investigated in connection with regenerative composition structures [HI [12] • 

Intuitively, a high concentration of the measure u near increases the intensity of 
small jumps of n„, hence increases the number of collisions X„. This is well exemplified 
by different kinds of the asymptotic behaviour of X„ for measures u with densities z/(dx) = 
cx°'~^{l — xY~^dx, with a,b,c> 0. For < a < 1 and b > the order of growth of X„ is 
linear [ISl [IS] , for a = 1 and 6=1 the order is n/ log ra [HI [3 [TS], |2I] , for 1 < a < 2 and 
6 = 1 the order is n^~° [12] , while for a = 2 and 6 > the order is log^ n [IB] . The class 
of beta distributions covered by the present paper corresponds to the range a > 2, when 
Xn exhibits a logarithmic growth. See [2] for other aspects of the beta-coalescents. 

We shall assume throughout that the characteristic measure z/ is a probability measure 
on the open interval (0, 1) such that 



(i) the support of v is not contained in a sequence of the form (1 — ^7", n = 0, 1, . . .) 
for some (5 > and 7 G (0, 1), 

(ii) V satisfies 

/ I logx|z/(da;) < 00. (2) 

i(o,i) 

Note that the more general case of a finite measure v can be reduced to the case of a 
probabihty measure by a hnear time change of the coalescent . 

2 Annihilator and a coupling 

Adopting Pitman [22l p. 1884], the coalescent process n„ can be constructed as follows. 
Let < Si < S2 < • • • be the arrival times of a unit Poisson process and let 1— r/i, 1— ''725 • • • 
be an iid sample from i/, also independent of the Poisson process. At time Si each of 
the n particles is marked 'head' with probability 1 — r^i or 'tail' with probability r^i. 
Those particles marked 'head' are merged in one particle, provided there are at least two 
'heads'. Thus k > 2 'heads' imply that n — k + 1 particles remain after merging, i.e. 
nn(si) = n — k + 1, and A; = 1 or = 'heads' imply that n„(si) = n. The procedure is 
then iterated at time S2 with the remaining particles and the new probability for 'heads' 
being 772, and so on. It is immediately checked that this description agrees with the rates 
([T]). The rules define the n„'s consistently for all values of n, a property used by Pitman 
[22j to introduce a partition-valued coalescent process with infinitely many particles. 

There is a simpler process on {0} U [n], whose transition rate from m to m — k for 
k = 1, . . . ,m is [!^\n,k (where A^j is defined by the same integral ([1]) with k = 1). Let 
K„ = (K„(t),t > 0) denote the process starting at n and terminating at 0. We shall call 
this process annihilator in view of the following interpretation: when there are m particles, 
any /c-tuple of them collides and annihilates at rate Xm,k, so that after the collision m — k 
particles remain (we stress that the case fc = 1 of 'singular' collision is not excluded). A 
sequential construction of the annihilator is almost the same as that of the coalescent. At 
the generic time Sj of a unit Poisson process a random value 1 — rjj is sampled from i/, 
and each of the remaining particles (if any) is marked 'head' with probability 1 — rjj or 
'tail' with probability t^j, then the particles marked 'heads' are removed. 

The sequence of positive decrements of K„ is a random composition (ordered partition) 
of the integer n. For instance, the path 10, 7, 2, 1, 0, ... of Kio corresponds to the composi- 
tion (3, 5, 1, 1) of the integer 10. As n varies, the random compositions are consistent with 
respect to a sampling operation, i.e. they comprise a composition structure, see P^ [TU] . 
Two features differ the transition rules of n„ and K„: 

(i) if at time Sj exactly fc > 2 of the remaining particles are marked 'heads' then n„ 
decrements by /c — 1, while K„ decrements by k, 

(ii) if A; = 1 then n„ has no transition at Sj, while K„ decrements by 1. 

Note that neither n„ nor K„ jump at time Sj if none of the remaining particles is marked 
'head'. 



To couple the coalescent and the annihilator we just define them using the same unit 
Poisson process with arrival times Si, S2, • • • and the same sample 1 — 771, 1 — 772, .. . from 
V. We shall call the initial n particles primary and their followers resulting from mergers 
secondary. The primary particles are labelled 1,2, ... ,n. If at some time exactly A; > 1 
particles are marked 'head', these k particles are replaced by one secondary particle. The 
variable K„(t) is identified with the number of primary particles alive at time t, and n„(t) 
is identified with the total number of primary and secondary particles at time t (thus by 
this coupling we always have n„(t) > K„(t)). For instance, if at time si the particles 
2,3 and 6 are marked 'head', the number of primary particles decreases by 3, and one 
secondary particle is born. 

Let Kn be the number of transitions of K„ as the process proceeds from n to 0, and 
let Kn,i be the number of jumps of size 1 (singular collisions). Introduce the absorption 
time of the annihilator (T„ = inf{t > : K„(t) = 0}. We define further Knfl to be the 
number of epochs Sj < (t„ when none of the primary particles is marked 'head'. We 
define Un to be the number of particles that remain at time o"„; all these particles are 
secondary, so f/„ = n„(o"„). For two real- valued random variables X and F, we say that 
X is stochastically smaller than Y and write X <(iY /li P{X < x) > P{Y < x) for every 
X eR. 

Lemma 2.1. The following stochastic order relations are true 

Kn -Kn,l <d Xn <d Kn + Knfi + Xu^, (3) 

where in the right hand side {Kn, Knfl, Un) is assumed to he independent of {Xj : j G N}. 

Proof. We use the above coupling. By collision of at least two primary particles both 
processes n„ and Kn decay, hence Kn — Kn^i < Xn. The number of jumps of n„ before 
and including time o"„ does not exceed the number of Poisson epochs, which is Kn + Kn^, 
and after time (T„ the coalescent evolves with f/„ particles, independently of the history 
(given Un)- This yields the upper bound for X„. D 

Lemma 2.2. The sequence of distributions of Xu„ {n G N) is tight. 

Proof. First of all, 1 < Xjj^ < Un because the number of collisions after o"„ does not 
exceed the number of particles that remain. To bound f/„, we focus on the behaviour of 
the secondary particles. At any time Sj < an, each of the secondary particles is removed 
with probability 1 — rjj and at most one secondary particle is born (possibly by collision 
of only primary particles). We can only increase the number of secondary particles if we 
assume that a secondary particle is always born at every time Sj < an. 

For j G No := N|J{0} denote by Qj the number of secondary particles at time Sj 
(so •= 0) and let JF be the o"-field generated by all the variables rjk, /c G N. Note that 
Qq = 0. With the above convention, given JF, the sequence {Qj : j E Nq} becomes 
a time-inhomogeneous Markov chain on the set of non-negative integers with transition 
probabilities 

P{Q, = A; I Q,_i = m} = f"^ X^-\l - r^^)^->'+\ k e {1, . . . ,m + 1}. 



Set 

Nn := inf {A; > 1 : nr]i • • • r^^ < 1}, n G N. 

Almost surely on the event {Nn > 2}, 

E {Qn„ I J^) = 1 + ^ r]jTlj+i ■ ■ ■ r/^„ =: 1 + e„, nE N. 
i=2 

Even though the subsequent argument is the same as in the proof of [20l Theorem 1.1], 
we include it for completeness. For fc = 2, 3, . . ., we have 

Onl{7V„=fc} = ?/l • • • ^iV„ I — + + . . . + I l{7V„=fe} < 

< \ h . . . H 1{JV„=A:} a.s., 

n V?! ^1^2 r]i...r]k-ij 



which implies that 

By the key renewal theorem 



1 ^^^ 1 



n 



:^ ^i---r/fc 



lim -E > = 

/"logn / °Q 

lim / e-(i°§"-^) y P{- logr/i - ... - logr^fe G dx} 



E(— logr]) ' 



vfc=0 

no matter whether mi := E(— log?]) is finite or infinite. Since EQAf„l{Af„=i} = o(l), we 
have proved that 

limsupE^AT^ < l/mi + 1, 

n— »oo 

which implies that the sequence of distributions of Q^^ {n G N) is tight. Note that A^„ 
is the number of steps that the multiplicative renewal process {jji ■ ■ -rij, j = 0, 1, . . .), 
needs to drop below level 1/n. Now we wish to replace the fixed level 1/n by a random 
level associated with the last primary particle disappearing at time an] this random level 
can be identified with the minimal order statistic i?i_„ of the sample of size n from the 
uniform [0, 1] distribution (see the remark after Theorem 13.11 for this connection). From 
the subadditivity of the first-passage time processes for random walks, and from the 
convergence in distribution of — log -Ri^„— log n (which is a standard fact from the extreme- 
value theory), it follows that the number of steps p„, say, the renewal process needs 
to cross the interval (min(i?i „, l/n), max(i?i,fi, 1/n)) remains stochastically bounded as 
n — >• cxD. On the event An := {Ri,n < 1/"^}) as the process proceeds from 1/n to Ri^n, the 
number of secondary particles may increase at most by pn- On the complementary event, 
Un — QN„ < XnUn a.s., whcrc Xn ^ (0, 1) is a stochastic sequence (depending on p„) which 
is bounded away from one. The latter inequality follows from the observation that in 
the coalescent process under consideration the number of particles involved in the generic 
collision is proportional to the number of existing particles. In conclusion, {Qnu ~ Un\ is 
stochastically bounded on both events, which completes the proof. D 



3 Number of collisions 

With these prehminaries, we are ready to derive the asymptotic behaviour of the number 
of colhsions X^ from the previously known results on Kn [12j. Introduce the logarithmic 
moments 

mi := E(- logr]), m2 := Var(log?7) 

for a random variable 1 — rj with distribution u and recall the standing assumptions of 
the paper given at the end of Section [TJ It is worth mentioning explicitly that in parts 
(l)-(3) of the next result mi < oo. 

Theorem 3.1. The following two assertions are equivalent: 

(i) There exist constants {an, &„ : ri G N} with a„ > and 6„ G M such that, as n -^ cxd, 
{Xn — bn)/cLn convcrgcs weakly to some non- degenerate and proper distribution. 

(ii) The distribution of {—logrj) either belongs to the domain of attraction of a stable 
law, or the function x t—^ P(— logrj > x) slowly varies at oo. 

Furthermore, there are five different regimes of convergence. 

(1) //m2 < oo, then, with bn '■= m^f^logn and an := (m|f'^m2 logn)"*^'^, the limiting law of 
{Xn — bn)/an is Standard normal. 

(2) //m2 = oo and 

(logy)^ z/(dy) ~ L(— logo;) as X ^ 0, 



(x,l) 



'gn\, 



for some L slowly varying at oo, then, with bn := m]""'^logn and an := m^ C[io, 
where Cn is any sequence satisfying limn^oo nL{cn) / c^ = 1, the limiting law of {Xn — 
bn)/0'n is standard normal. 

(3) Assume that for some function L slowly varying at oo, the relation 

< x) ~ (— logx)~°L(— logx) as X — s> 0, (4) 



holds with a E [1,2), and that rai < oo if a = I. Then, with bn '■= m^^^logri and 



m 



-{a+l)/a 



C[iogn]; where c„ is any sequence satisfying \im.n^oo nL{cn)/c'^ 



the limiting law of {Xn — bn)/an is a-stable with characteristic function 

t ^ exp{-|t|°r(l-a)(cos(7ra/2) + zsin(7ra/2)sgn(t))}, t G M. 

(4) Assume that mi = oo and that (^ holds with a = 1. Let c be any positive function 
satisfying lim^^^xL{c{x))/c{x) = 1 and set ^{x) := x J^^^^_^^^y^¥{r] < y)/ydy. 
Let b{x) be any positive function satisfying b{ip{x)) ~ ip{b{x)) ~ x. Then, with 
bn '■= b{\ogn) and an '■= b(\ogn)c{b{\ogn)) / \ogn, the limiting law of {Xn — bn)/an 
is 1-stable with characteristic function 

t ^ exp{-|t|(7r/2-2log|t|sgn(t))}, teR. 



(5) // (jl]) holds with a G [0,1), then, with a„ := log"n/L(log?7,), the limiting law of 
Xn/dn is scaled Mittag-Leffter 6a (exponential, if a = 0) with moments 

x^'Oaidx) = —, -— -, nGN. 

[o,oo) r"(l -a)r(l + na) 

Proof. Following the same line of argument as in the proof of Theorem 2.1 in p2|, one 
can show that a necessary condition for the convergence of [X^ — bn)/an is that the 
scahng satisfies a„ — »• cxo, which will be assumed in the sequel. To derive the result 
from Theorem 2.3 in [12] it is enough to show that (X„ — &„)/a„ weakly converges to a 
proper and non-degenerate distribution if and only if {Kn — bn)/an weakly converges to 
the same distribution. This equivalence is implied by Q. Indeed, by Lemma [2T2| Xjj^/an 
converges in probability to 0. Under the assumption ([2]), Kn^ converges in distribution 
by [121 Theorem 2.2], and Kn^ converges in distribution by [121 Proposition 5.2] (see also 
[TT] for the case mi < oo). D 

Remark. We describe briefly a more general construction with possibly infinite charac- 
teristic measure z/ satisfying j,^ ^. xuldx) < oo. Let u be the image of u under the transfor- 
mation X H-^ — logx, let (St) be a subordinator with Levy measure u, and St = exp{—St). 
Then {St) is a non-increasing positive pure-jump process with 5*0 = 1. Primary parti- 
cles are labelled by independent uniform[0, 1] points i?i, . . . , i?„, also independent of (St). 
The annihilator process K„ = (K„(t),t > 0) is defined by setting K„(0) = n and by 
decrementing the process at time t by fc if the interval [St, St-] contains exactly k of the 
uniform points. Thus K„ may have a jump at t only if St- > St- The coalescent n„ 
is defined in a similar way, but except Ri, . . . ,Rn representing primary particles, new 
points representing secondary particles appear. Set n„(0) = n. Conditionally given that 
[St, St-] contains k > 1 points, primary and secondary, let Un decrement by A; — 1 and let 
one point with uniform distribution on [0, St] appear to represent a secondary particle. 
The inequality Xn > Kn — Kn,i is still true, but another bound in ([3]) makes no sense for 
infinite u, since a variable analogous to Knfl is identically oo. It seems plausible that in 
this more general setting X„ and Kn — i^n.i still have the same limit laws and require the 
same scaling. Asymptotics of Kn, Kn^i were studied for infinite measures z/ satisfying a 
condition of regular variation [T3l |1] 



Example 3.2. Application of case (1) of Theorem 13.11 to coalescents governed by beta 
measures u^dx) = cx'^~^(l — x)^~^dx with parameters a > 2 and b,c > shows that the 
scaled and centered number of collisions 



(X„ - m^ ^ log n)/{ni2m-^ ^ log 



n 



,1/2 



converges in distribution to the standard normal law, where mi = \E'(a — 2 + 6) — \t'(6) and 
m2 = \1''(&) — \&'(a — 2-1-6), and \E' is the logarithmic derivative of the gamma function. 

For the special case uldx) = dx (uniform distribution) we have gn,m/gn = l/(^ ~ !)• 
Thus, n„ jumps from n to every m & {1, . . . ,n — 1} with the same probability. In this 
case Theorem 13.11 implies that (X„ — log n) / \/log n is asymptotically normal, which agrees 
with a classical result known in the context of records or cyclic partition of permutations. 
We refer to the survey of [1] for some further information along with references. 



More generally, when z/ is a beta(l, b) distribution, a straightforward calculation shows 
that n„ jumps from n > 2 to any state m E {1, ... ,n — 1} with probability Wm/{wi + 
• • • + Wn-i), where Wk := T{k + b — 1)/T{k). By p5^j Proposition 2], X„ is representable 
as the sum of n — 1 independent indicators Ik with 

P(4 = l) = ^ "^^ , ke{l,...,n-l}. 

It can be shown that P(/a; = 1) ~ b/k as k ^ cx). From this, as n — >■ cxd, the expectation 
and the variance of X„ are both asymptotically equal to 61ogn. By the central limit 
theorem, the distribution of (X„ — 6 log ra) /a/6 log n converges to the standard normal 
law, which again agrees with Theorem 13.11 case (1). The latter result is classics which 
dates back to [T7] . 

The merging evolution of n particles in the coalescent is often represented by a rooted 
tree with n leaves, and with the internal nodes corresponding to collisions [3], |23]. One 
interpretation of the similarity of the asymptotic behaviour of X„ and Kn is that almost 
all internal nodes of the tree are directly connected to one of the leaves. 

Some other properties of the coalescent follow easily in the same way. Let Z„ be the 
time of the first collision of a tagged particle (say, labelled 1). In the tree representation, 
Zn is the length of an external branch, which connects a leaf to its closest internal node. 
The asymptotic behaviour of Zn for the Kingman coalescent and the Bolthausen-Sznitman 
coalescent (z/(dx) = x~^dx) was studied in [5] and [8], respectively. In our setting we have 
the following. 

Proposition 3.3. As n ^ oo, Zn converges weakly to an exponential variable with pa- 
rameter p := E(l — ry) = J,g ^. xuldx) . Moreover, the number of collisions in the coalescent 
before and including time Zn converges weakly to a geometric variable with parameter p. 

Proof. In the process with infinitely many particles a tagged particle is marked 'head' 
at rate p, and when this occurs the particle collides with infinitely many other particles, 
because the assumption that z/ is a probability measure implies that the coalescent does 
not come down from infinity [22j (which is also obvious from the construction). This 
proves the first claim. The second claim follows similarly, by noting that the time of the 
first collision of a given primary particle is Si with probability p. D 

4 Absorption times 

The coupling allows us to derive the asymptotics of r„. 
Lemma 4.1. The following stochastic order relations hold 

0"n <d Tn <d Cfn + t'jj^-, 

where, on the right-hand side, t'- = Tj, j G N, and an, (tJ) and Un are assumed to be 
independent. 



Proof. These relations are obvious from the couphng. When all the primary particles 
disappear at time an, there are Un secondary particles left. D 

From the results obtained in [12] it follows that an has the same limiting distribution as 
the first passage time through the level logn for a compound Poisson process {Rt '■ t > 0} 
with the generic jump (— logr/) and intensity one. By Lemma [2.21 and Lemma [4.11 the 
same asymptotics is valid for ?:„. This argument leads to the following result. 

Theorem 4.2. Assume that condition (ii) of Theorem \3.1\ holds. Then, with X„ replaced 
by Tn, condition (i) and parts (2) - (5) of that theorem hold without changes, while part 
(1) holds with an = (mf^(m2 + mf) logn)-*^/^. 

Proof. We focus on the asymptotics of the first passage time process {Ux '■ x > 0} defined 
as 

Ux := inf{t >0:Rt>x}, x>0. 

Let Nx be the number of renewals on [0, x] of the renewal process with A^q = and steps 
distributed as (—log?]). Obviously, Ux = sn^, where as above si,S2,--- are the arrival 
times of an independent unit Poisson process. Assume that condition (ii) of Theorem 13.11 
holds, and let a{x) and b{x) be the same as in p^J Proposition 8.1] (setting a„ := a(logn) 
and bn '■= b{logn) gives the normalizing constants for the cases (2) - (5) of Theorem 13. II) . 
We will exploit the following decomposition 

Sat, - b{x) sn^ - Nx y/Nx , Nx - b{x) 



a{x) \/Nx cl(x) a{x) 

In view of the independence and since A^^. —^ cxd almost surely as x — ;> cxo, we have by the 
central limit theorem that, as x ^ cx), {sn^ — Nx)/\/7^ converges in distribution to a 
random variable ^i with the standard normal distribution. By |12. Proposition 8.1], the 
right-most term in ([5]) converges in distribution to a random variable ^2 with a proper 
and non-degenerate law. Note that ^1 and ^2 are independent. 

Case (1). By [121 Proposition 8.1], b{x) = m^^x, a{x) = (m2mj"^a;)^/^, and ^2 has the 
standard normal law. Since, by the strong law of large numbers, as x — > 00, Nx/x — ^ m]"^ 
almost surely, we conclude that the left-hand side in ([5]) converges in distribution to a 

— 1/2 

random variable mim2 ^1 + ^2 which has the normal law with zero mean and variance 
m^m^"^ + 1. Therefore, (Ux — mj"^x)/(ni]^'^(m2 + m^)x)^/^ weakly converges to the standard 
normal law, as desired. 

Cases (2) - (5). By [T?, Proposition 8.1], as x — > cx), \fNl/a{x) converges to in 
probability. Therefore, the left-hand side in ([5]) weakly converges to the same law as the 
right-most term in ([5]). We again refer to [121 Proposition 8.1] for the identification of the 
normalizing constants and the relevant limiting laws. D 

Example 4.3. Theorem 3.2 is obviously applicable to the class of beta-coalescents with 
parameter a > 2 and 6, c > considered in Example 13. 2[ Consequently, for such beta- 
coalescents, the scaled and centered absorption time (r„— mj^^ logr2)/(nij"^(m2-|-m^) logn)^/^ 
converges in distribution to the standard normal law, with mi and m2 as given in Example 
[33 
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